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An approach to general relativity based on conformal flatness and quasiequilibrium (CFQE) as- 
sumptions has played an important role in the study of the inspiral dynamics and in providing 
initial data for fully general relativistic numerical simulations of coalescing compact binaries. How- 
ever, the regime of validity of the approach has never been established. To this end, we develop an 
analysis that determines the violation of the CFQE approximation in the evolution of the binary 
described by the full Einstein theory. With this analysis, we show that the CFQE assumption is 
significantly violated even at relatively large orbital separations in the case of corotational neutron 
star binaries. We also demonstrate that the innermost stable circular orbit (ISCO) determined in 
the CFQE approach for corotating neutron star binaries may have no astrophysical significance. 
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a. Introduction The analysis of general relativistic 
binary neutron star processes is an important yet chal- 
lenging endeavor. The importance of studying binary 
neutron star coalescence is rooted in observational as- 
tronomy; the process is considered to be a candidate 
for the sources of detectable gravitational waves as well 
as X- and 7-rays. Recently, there has been a large ef- 
fort expended towards numerically simulating the inspi- 
ral and coalescence of compact objects in order to deter- 
mine their signature gravitational waveforms in prepara- 
tion for the detection of gravitational waves by the new 
generation of gravitational wave observatories, including 
LIGO, LISA, VIRGO, TAMA, and GEO. The complexity 
of the nonlinear Einstein field equations, the lack of ex- 
act symmetry of the system, and the coupling of different 
timescales in the process combine to make the study of 
the neutron star inspiral coalescence a major challenge. 

A recent approach used to investigate the inspiral coa- 
lescence problem that has drawn much attention P, [E 
IE EL IE IE IE III Ell is based on the separation of 
two timescales, namely, the orbital motion timescale and 
the gravitational radiation timescale. This approach, 
which we refer to as the conformally flat quasiequilib- 
rium (CFQE) approach, seeks to construct general rela- 
tivistic configurations corresponding to two compact ob- 
jects in circular orbit by assuming that the spacetime is 
essentially stationary with an approximate Killing vec- 
tor (quasiequilibrium, QE). The slow secular evolution is 
driven by the gravitational radiation which has a much 
longer timescale than the orbital motion timescale, so 
that the spacetime is thought to be describable by a con- 
formally flat (CF) spatial metric. 

The CFQE approach leads to a two-phase description 
of the inspiral and coalescence process. In the "slowly 
inspiraling phase", the evolution is approximated by a 
time sequence (which we refer to as a CFQE-sequence) 
of CFQE configurations with the same rest mass and spin 
state (e.g. corotating or nonrotating). Each CFQE con- 
figuration in the time sequence is obtained by solving the 



four constraint equations plus one of the dynamical Ein- 
stein equations under the CF and QE assumptions (the 
QE assumption affects the constraint equations only indi- 
rectly) . The inspiral is assumed to proceed along the se- 
quence in a quasistationary manner on a timescale much 
longer than the orbital timescale. The slow secular evo- 
lution is driven by gravitational radiation, which can be 
calculated using e.g., the standard quadrupole formula 
(see, e.g., ^3)- ^ n the corotating case, when the orbit 
shrinks to a certain radius (which will depend on the 
equation of state (EOS) of the neutron stars), the CFQE 
configuration becomes secularly unstable, with the ADM 
mass of the system attaining a minimum among all con- 
figurations in the sequence. The orbit at this radius is 
referred to as the innermost stable circular orbit (ISCO). 
One expects the binary to become dynamically unstable 
at an orbital radius that is slightly less than this ISCO 
point. This ISCO point is therefore thought to be the 
point at which the system enters the "plunge phase" , in 
which the subsequent merging occurs within an orbital 
timescale. 

Besides providing a description of the inspiral coales- 
cence process, the CFQE approach is important for an- 
other reason: it is expected to provide a starting point 
for numerical investigations in general relativity. The 
setting of initial data is a sensitive issue in the numer- 
ical integration of the Einstein equations. On the one 
hand, astrophysically relevant initial data must be used 
in order to make contact with gravitational wave obser- 
vations. On the other hand, we are forced to use ini- 
tial data in the late state of the inspiral, due to limited 
computational resources. As the CFQE approach gives 
configurations that satisfy the constraint equations of the 
Einstein theory, one is tempted to start numerical inte- 
gration with a CFQE configuration near the ISCO point 
as determined in the CFQE approach, see, e.g., ^^|- The 
danger in using the CFQE approach to determine the 
initial configuration for numerical integration is the un- 
known astrophysical relevance of the initial data, regard- 
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less of the fact that each CFQE configuration satisfies 
the constraints of general relativity and is therefore, in 
principle, a legitimate initial data set for numerical inte- 
grations. In order for the results of the numerical evolu- 
tions to be relevant to observations, e.g., the gravitational 
waves emitted in neutron star coalescences, we have to 
make sure that the initial data actually corresponds to a 
configuration in a realistic inspiral. The setting of physi- 
cal initial data is a standard difficulty in numerical rela- 
tivity and it is particularly acute in this case; if the initial 
data is picked to be a CFQE configuration in a regime 
where the CFQE-sequence is not expected to be valid, 
the initial data must be regarded as highly contrived and 
the resulting calculation will have no astrophysical sig- 
nificance. 

In this paper we analyze the limitations of the CFQE 
approach, which suffers from two basic problems. The 
first problem is that the CF and QE assumptions violate 
the Einstein equations which, for this system, forbid both 
an exact Killing vector and a conformally Hat metric at 
all times (the metric can only be chosen to be confor- 
mally flat at one point in time but not along the entire 
time sequence). The second problem is that the CFQE 
approach itself does not provide a way to estimate the 
error involved in its assumptions. A determination of its 
accuracy has to be carried out with a completely differ- 
ent analysis. Because of this, the regime of validity of 
the CFQE approach has never been determined. Here, 
we calculate the rate of increase of the accuracy of the 
CF and QE approximations with increasing initial neu- 
tron star separation. For the corotating binaries studied 
(with a r = 2 polytropic equation of state), we demon- 
strate that an initial separation of at least 6 neutron star 
radii is needed in order for the violation of the QE as- 
sumption to remain below 10% (in some reasonable mea- 
sure defined below) after a short evolution (a fraction of 
an orbit). (The violation of the CF assumption is an or- 
der of magnitude smaller at this separation.) It is safe to 
conclude that, starting at a separation of 6 neutron star 
radii, the evolution (e.g., the orbital phase) of the system 
is completely different from that of the CFQE-sequcncc 
well before the CFQE ISCO separation (less than 3 neu- 
tron star radii) is reached. 

In what follows, we present the results of our analysis 
of the CFQE-sequence approximation for binary, coro- 
tating neutron stars. The details of all calculations and 
the subsequent analysis of the numerical results, along 
with results on long timescale (e.g., multiple orbital pe- 
riod timescale) numerical evolutions will be reported else- 
where 1 1 3| - 

b. A General Relativistic Analysis of the CFQE 
Approximations for Corotating Neutron Star Binaries. 
To analyze the validity of the CFQE-sequence approxi- 
mation, we compare it to solutions of the full Einstein 
equations whose initial Cauchy slice corresponds to spe- 
cific CFQE configurations. We construct these corotat- 
ing CFQE configurations following the algorithm detailed 
in [lj . The neutron stars we use are described by a poly- 



tropic EOS, P = Kp r , with r = 2 and K = 0.0445c 2 /p n , 
where p n is nuclear density (2.3 x 10 14 g/cm 3 ). For this 
EOS, the maximum stable static neutron star configura- 
tion has an ADM mass of 1.79M Q and a baryonic mass 
of 1.97Mq (Mq is 1 solar mass). In this paper, we use 
neutron stars whose baryonic mass Mo is 1.49M Q , which 
is approximately 75% that of the maximum stable con- 
figuration. The ADM mass of a single static neutron star 
for this configuration is 1.4M Q . 

The corotational CFQE configurations with a fixed 
baryonic mass can be uniquely specified by the separa- 
tion between the two neutron stars. We measure the 
spatial separation on a constant time hypersurface by 
the spatial geodesic distance £1.2 between the points of 
maximum baryonic mass density in each of the neutron 
stars. This particular choice of measure of distance de- 
pends on the choice of time slicing. In order to compare 
with the CFQE-sequence approximation in an invariant 
manner, we use the same slicing condition (maximal slic- 
ing) as that used in CFQE-sequence approach. Four 
CFQE configurations denoted respectively NS-1 to NS-4 
are used as initial data in our general relativistic simula- 
tions in this paper: £i t2 /M = 23.44,25.94,29.78,35.72. 
The corresponding orbital angular velocity parameters 
are QM = 0.01547,0.01296,0.01022,0.00746. 

In the following we first take care to construct invariant 
measures that can be used to compare, and give a sense 
of the difference between, two different spacetimes, one 
of which is obtained by solving the Einstein equations 
and the other represented by the CFQE-sequence, which 
does not satisfy the Einstein equations. 

c. The QE approximation. One basic assumption of 
the CFQE approach is the existence of a timelike, heli- 
cal Killing vector field. In the case of a corotating bi- 
nary, the 4-velocity u a of the fluid is proportional to 
this Killing vector. This implies the vanishing of the 
type-(0,2) 4-tensor Q ab = V ' a u b + V b u a + u a a b + u b a ai 
where a a = u h V b u a is the 4-acceleration of the fluid 
(V a denotes the covariant derivative operator compat- 
ible with the 4-metric g a b)- The quantity Q a bU a van- 
ishes identically. We can thus monitor the space-space 
components of Q ab during full general relativistic simu- 
lations as a way of monitoring how well the 4-velocity 
u a stays proportional to a Killing vector field. Define 
Qij as the projection of Q ab onto the constant t spa- 

tial slice: Q tJ = P l a P j b Q ab = Q ab (-£z)° '(^j)' where 
Pab = 9ab + n a n b is the projection operator onto the 
constant t spatial slices with unit normal n a . The 3- 
coordinate independent norm \Qij\ of this matrix Qij is 
the square root of the largest eigenvalue of Qij Q 3 k , where 
we have raised and lowered 3-indices with the 3-metric. 

To gain a sense of the size of \Qij\, we normalize 
\Qij\ by the norms of its two principle parts: Q = 
|Q i j|/max{|<9i ij | , |<92ij|}, where Qi ab = V a u b + W b u a , 
and Q-2 a b = u a<^b + u b a a . A value of Q = signifies that 
the 4-velocity of the fluid is proportional to a timelike 
Killing vector as in the CFQE-sequence, while a value of 
Q = 1 would signify that the assumption is maximally 
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FIG. 1: The quantity (Q) is plotted as a function of time 
for fully consistent general relativistic numerical calculations, 
using CFQE configuration NS-1 to NS-4 as initial data. 
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violated. 

Since Q is meaningful only inside the fluid bodies, a 
natural global measure of the magnitude of Q is its the 
baryonic mass weighted integral, denoted by (Q): 

where the integrals are taken to be over the entire spatial 
slice, with support only where p / 0. Here, we have 
denoted 7 to be the determinant of the 3-metric and W 
to be the Lorentz factor of the fluid so that J d 3 x y/jpW 
is the total (conserved) rest mass of the system. 

In Fig. we plot (Q) as a function of time in the 
general relativistic numerical simulations, using CFQE 
configurations NS-1 to NS-4 as initial data. For each ini- 
tial CFQE configuration, the subsequent solution to the 
Einstein equations has (Q) reaching a maximum value 
after a very short time (between 0.5% and 1.5% of an or- 
bit, 27r/f2 being the orbital period corresponding to each 
CFQE configuration). These results arc in contrast to 
the CFQE-sequence approximation, where (Q) is exactly 
zero. The rapid growth of (Q) to such a large value sig- 
nals that, in the present case of corotating neutron stars, 
after just a short time the CFQE-sequence is very dif- 
ferent from the spacetime described by the full Einstein 
equations even when the latter is started with the same 
initial data. Within a small fraction of an orbit, the 4- 
velocity of the fluid already is in no way proportional to 
a Killing vector. 

Also, we see from Fig. ^ that the maximum value ob- 
tained by (Q) decreases from approximately (Q) = 0.27 
to (Q) = 0.17 as the geodesic separation of the initial 
CFQE configuration increases from ^1/2/A^o = 23.44 to 
Ix/i/Mo = 35.72. In principle one can obtain an arbi- 
trarily accurate CFQE configuration by starting the in- 
spiral evolution with a large enough separation. While 
starting at a larger separation does not imply the CFQE- 
sequence will be more astrophysically relevant at late 
times, the fully relativistic simulation starting with it 



FIG. 2: Fig. 2a shows the maximum value of (Q) obtained 
in our fully consistent general relativistic simulations using 
CFQE configurations with different geodesic separation £1,2 
as initial data. Fig. 2b shows the corresponding maximum 
value of [H ) 



will be. This effect can be seen directly in Fig. 2a, where 
the maximum value of (Q) obtained in the general rela- 
tivistic simulations is plotted as a function of the initial 
geodesic separation t\^. A careful analysis of the numer- 
ical errors (shown as error bars in the figure) based on a 
Richardson extrapolation technique (for both truncation 
errors and boundary errors) can be found in In this 
way, we can determine the required separation of the ini- 
tial CFQE configuration such that the QE assumption 
is valid to any required level in our fully consistent gen- 
eral relativistic simulations. In Fig. 2a, we fit a reciprocal 
power law, jp^- + ,, a2 . y + rr^s + ic ai ^ , to the maximum 

l«l,2j (1-1,2) ('1,2) (A j2 ) 

value attained in our general relativistic calculations (as 
well as to the lower bound of the estimated error). We 
can see that, e.g., one would have to use a CFQE configu- 
ration initial data set with a geodesic separation between 
the neutron stars of about £1$ = 47 Mq in order for the 
subsequent solution to have the fluid 4- velocity following 
a Killing vector to the 10% level. 

d. The CF Approximation. While the 3-metric can 
be chosen to be conformally flat on any one time slice, the 
assumption that the 3-metric remain conformally flat for 
all time violates the dynamical Einstein equations. Again 
we begin with an invariant construction of a measure of 
the violation. Because conformal flatness is a property 
of a constant time slice, and because we are using the 
same slicing conditions in our simulation as that used in 
the CFQE-sequence approximation, we seek a 3- invariant 
(with respect to spatial coordinate transformations) to 
monitor the CF assumption. The Bach 3-tensor Bijk = 
2D^ (( 3 ^Rj] k — jjj^^Ft) can be shown to vanish if and 
only if the 3-metric 7^ is conformally flat. The Cotton- 
York tensor, Hij , is related to the 3-Bach tensor by Hfj = 
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FIG. 3: The quantity {H) is plotted as a function of time 
for fully consistent general relativistic numerical calculations, 
using CFQE configuration NS-l to NS-4 as initial data. 



e mn jB mni , where e«fc is the natural volume element 3- 
form. We define the scalar H as the matrix norm of Hij , 
normalized by the size of the covariant derivative of the 
3-Ricci tensor: H = \ H mn \ / ^JV, ^R 3k V l ( 3 )i# fe , where 
\Hij \ denotes the matrix norm of the As before, we 
define the baryonic mass density weight of H, denoted as 



(H) = 



J d 3 x \H\ ^pW 
Jd^x^pW ' 



(2) 



where the integrals are taken to be over the entire spatial 
slice, but have support only where p 7^ 0. 

In Fig. |3 we plot (H) as a function of time for the 
general relativistic simulations using initial data CFQE 
configurations NS-l through NS-4. The profile of the vio- 
lation of the conformal flatness assumption (H) is similar 
to that of the Killing vector assumption (Q) : a quick rise 
to a maximum value at a timescale on the same order of 
but slightly longer than that of (Q). This is contrasted 
to the CFQE-sequence approximation, which has (H) ex- 
actly equal to 0. 

In Fig. 2b, we plot the maximum value of the quan- 
tity (H) function of initial geodesic separation 



attained in our general relativistic numerical simulations 
using the four CFQE configurations NS-l through NS-4 
as initial data (see [b| for a detailed analysis of the trun- 
cation and boundary errors in our calculation, which we 
indicate as error bars in the figure). As in Fig. 2a, we 
fit a reciprocal power law to both the maximum value 
of (H) and the lower bound of the error. We see that 
a CFQE configuration with geodesic separation £1.2 of 
approximately 47 Mq or greater must be used as initial 
data in order for the full solution to the Einstein field 
equations to have a value of (H) to be 0.01 or less. 

e. Conclusions In this letter we present a method 
for analyzing the regime of validity of the CFQE-sequence 
approximation. We apply this method to the CFQE- 
sequence approximation of corotating binary neutron 
stars. We show that for geodesic separations less than 
roughly 45 Mo, which is about twice the geodesic sep- 
aration of the ISCO configuration, the QE approxima- 
tion is severely violated and the CFQE sequence cannot 
be taken as a reasonable approximation of the binary 
evolution for these separations. Numerical simulations 
starting with CFQE configurations as initial data with 
geodesic separations smaller than 45 M$ cannot, there- 
fore, be considered as approximating a realistic neutron 
star binary inspiral. 
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